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Abstract 

We analytically calculate to second order the correction to the asymptotic form of quasinormal 
frequencies of four dimensional Schwarzschild black holes based on the monodromy analysis pro- 
posed by Motl and Neitzke. Our results are in good agreement with those obtained from numerical 
calculation. 
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1 Introduction 



Quasinormal modes (QNMs) were originally observed in considering the scattering or emission of 
gravitational waves by Schwarzschild black holes pQ. It was found that a characteristic damped 
oscillation, which only depends on the black hole mass, dominated the time evolution in a certain 
period of time. Since then QNMs have been investigated extensively both analytically and numerically. 
For a general review and classification, see Refs. [2, 3J. From numerical studies, an asymptotic formula 
for quasinormal frequencies of Schwarzschild black holes was obtained [I] : 



spectrum of quantum black holes proposed in Ref. [5] . This was confirmed later by Motl and Neitzke 
[7J. The recent surge of interest in the QNMs derived from its possible application in determining 
the Immirzi parameter in loop quantum gravity [5]. The numerical value In 3 in the real part of the 
asymptotic quasinormal frequencies in Schwarzschild black holes was at first taken as a hint that the 
relevant gauge group in loop quantum gravity is SO(3) instead of the commonly believed SU(2). 
However, as shown in Ref. [7J, the value In 3 is not universal and one should take the argument with 
a grain of salt. 

Another interesting application of QNMs was pointed out by Horowitz and Hubeny in their study 
of a scalar field in the background of a Schwarzschild anti-de Sitter black hole [9]. According to 
AdS/CFT correspondence, a large black hole in AdS spacetime corresponds to a thermal state in CFT 
[TO] , They argued the decay of the scalar field corresponds to the decay of a perturbation of this state. 
In the BTZ black hole, a one-to-one correspondence was found between the QNMs in the bulk and 
the poles of the retarded correlation function in the dual conformal field theory on the boundary |llj . 
The idea of dS/CFT correspondence has also been proposed and formulated [15]. Since there is a 
cosmological horizon in de Sitter spacetime, QNMs may also be defined in principle. Similar studies of 
QNMs have also been carried out in de Sitter spacetime trying to lent support for such correspondence 
[16]. However, the situation there is more subtle and it seems QNMs only exist in odd dimensions [3]- 
Therefore, it is not clear whether such correspondence makes sense in even dimensions, and further 
study is necessary. 

2 Pert ur bat ive calculation of the asymptotic form of quasi- 
normal frequencies 

In Ref. |12j . the author calculated the first order correction to the asymptotic form of quasinormal 
frequencies of a Schwarzschild black hole using a WKB analysis. The result was extended to include 
the scalar field case using the monodromy analysis developed by Motl and Neitzke [13] . The agreement 
with numerical results is excellent. We will begin with a brief review of their method which made 
systematic expansion more accessible. In a background spacetime described by a metric , a massless 
scalar $ satisfies the following Klein-Gordon equation: 




(1) 



The real part in the above formula was later postulated to be In 3 [5] based on a discrete area 




(2) 
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For four dimensional Schwarzschild black holes, the metric is given by 

ds 2 = -f(r)dt 2 + /(r^dr 2 + rW, 
with f(r) = (1 - ^) and r = 2GM. Let 

$(r,t,n)=r<p(r)Y lm (tye iut . 
<j){r) now satisfies the following equation: 

-/(»-)J: /W^ + v(*O0 = «>V> 

with 



(3) 



(4) 



^) = (1--) 
r 

By a simple modification in the potential V(r) [2J, 



*(* + !) | r 



V(r) = (1 - ^) 



1(1+1) , (1-J 2 )r 



(5) 



the previous equation can also describes linearized perturbation of the metric or an electromagnetic 
test fields. Here, j = 0,1,2 which is the spin of the relevant field. They can also be classified as 
the tensor, vector, and scalar types of perturbation to the background Schwarzschild metric using the 
master equations derived by Ishibashi and Kodama [13]. Introducing the tortoise coordinate: 

x(r) = r + r ln(r/r - 1), 

one obtain a Schrodinger-like equation 



dx 2 



+ V[r{x)] \<f> = uj 2 <\> 



(6) 



Because of our convention in eq ©, QNMs are defined through the following out-going wave 
boundary condition: 



assuming Reoj > 0. Define 

f{x) = ^ 



as x — > — oo (horizon), 
asx^oo (spatial infinity), 



1 



-oo, 
as x — > oo. 



(7) 



(8) 



According to Ref. [7], the boundary condition at the horizon translates to the monodromy of f(x) 
around it 



M(r ) 



a 47rwro 



(9) 



The same monodromy can also be accounted for by those around r — and r — oo, and it has 
been shown that only the former one is non-trivial. To find the monodromy around r — 0, one need 
to introduce the complex coordinate variable 



z = uj(x — inro) — io[r + tq ln(l — r/fo)], 



(10) 
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which is vanishing at the black hole singularity r = 0. In the limit |r/ro| <C 1, the potential can be 
expanded as a series in y/ z/(ujro): 

v(z)= ^ 2 (1-J 2 ) 3l(l + l) + l-f 3l(l + l) + l-f 
1 j 4z 2 6^(-wr )V2 z 3/2 36V2(-wro) 3 / 2 z 1 /2 1 j 

Note that the third term in the above expression is of order (— wro)~ 3 / 2 and would not contribute 
until we consider third order perturbation. To second order in perturbation theory, the wavefunction 
can be expanded as 

+ ^=^« + —4> m + o(^- 3/2 ). (12) 



The zeroth, first and second order equations are given by 





dz 2 


V 4z 2 






dz 2 + 


^ 4z 2 


d^( 2 ) 


fl-J 2 


<fe 2 + 


V 4z 2 



respectively. Here, 



Define qf>^ (z) to be the two linearly independent solutions to the zeroth order equation 



+ 1J (o) = 0; (13) 
+ l) </> (1) =V=^8V(z)(f ) W; (14) 
+ l) </> (2) = V=W)5V(z)(l> w , (15) 



<l>f(z) = J — J ±j/2 {z). (17) 



In the asymptotic region z» 1 

4 O) (20«COs[z-7T(l±i)/4]. (18) 

It has been shown by Musiri and Siopsis that <j)±^ can be expressed in terms of cp± 

$\z)=C<l>f(z) [ Z dz x 8V{z 1 )<l> { °\z 1 ) ( l%Xz 1 )-Cc^\z) f Z dz 1 SV( Zl )^ ) (z 1 )^ ) ( Zl ); (19) 



ct&{z)=C^\z) dz l 8V{zi)^ ) {z l )4> { °\z 1 )-C^\z) dzx8V{z l )<t>f{zi)<f> { °\z l ). (20) 
Jo Jo 

where C = \/—u)T§l sin(7rj'/2) [13] . Similarly, <j>^ can in turn be expressed in terms of <fi±^ 

$\z)=C$\z) f dz 2 8V{z 2 )^\z 2 )<j^\z 2 )-C4 { -\z) f dz 2 5V{z 2 )cj ) f{z 2 )^ ) {z 2 ); (21) 
Jo Jo 

(2) (z) = C<$ (z) [ Z dz 2 SV(z 2 ) 0L O) (z 2 ) W (z 2 ) - Ctf& (z) f dz 2 SV(z 2 ) (z 2 ) 4> W (z 2 ). (22) 
Jo Jo 

In the limit, z — > oo, 

4 1} (z) = c_± <^ 0) (z) - c +± 0<_°>(z); (23) 
4 2) (z) = d_ ± (*) - d +± ^ (z). (24) 
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Here, 



c ±± =C I dz l 8V{z 1 ) ( j> < i\z 1 ) l j>f{z 1 ); 



(25) 



d±± = C 2 



dz 2 I dz 1 SV(z 2 )SV(z 1 ) ( f>^(z 2 ) (*a)0L O) (*i) ~ ^ ) (z 2 )^ ) {z 1 )\ 4 0) (zi). (26) 



Notice that <f>± defined in eq (| 1 T[) arc in fact linearly dependent to each other when j is an even 
integer. As a result, each of these coefficients is divergent by itself in these cases. It is reassuring to 
see that all the divergent pieces cancel among themselves so that physically interested quantities do 
have a smooth limit when j is an even integer. In zeroth order, the combination 



in the asymptotic region z 1. This can be extended to second order 

<t>{z) = 



(27) 



*?>(*)+ 



1 «#>(*) + — 



-wr 



-i7r(j/2) 



1 - 



C 



1 



W (z) + 



-ujr 



<t£\z) 



(28) 



by introducing two parameters £ and £. Naturally, they are determined by the condition that the 
coefficient of the e lz term is vanishing when z —>■ oo: 



£ = £++£-; 

C = + d ++ e i7rj '/2 _ d+ _ + e -i^/2 _ 



where 



C-| , £_ = c__ e 



Substitute the above result back to eq (|28|). we have 

£_ , £(£ m + c + _)-d__e-^/ 2 +d_. 



(z) = ie 47r ( 1 -^/ 4 sin(7rj72)e- 



1 



(29) 
(30) 



(31) 



(32) 



where the identity c |_ = has been used to simplify the expression. 

When going around the black hole singularity by 3ir, cj)± and <j}± both pick up an extra phase: 



z) = e 3i < 2± ^/ 2 4> { l\-z); 
<f,V>(e 3 ™z) = e 3 "( 3± J')/ 2 ^ 2) (-z). 



(33) 
(34) 



Consequently, 

0(e 3 ^z) = e 3i < 1+j y 2 

_ e -in(j/2) 

„3i7T(l-j)/2 



bf{-z)-i^L=<j>f{-z) 



-ujr 



<j>f(-z) 



1 



< 



\/-wr -wr 



(35) 
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To second order, 



[e 3in z) = -ie l,r(1 - i)/4 sin(37rj72)e- 



1 t (l+»e 3 ''J)g + + (l + »)g_ 

-(1 + i)££_ + [(1 + e i37rj )(d ++ e M ^ 2 - + 2d__ e -^/ 2 - 2d+_] 
+ -o;r (-l + c' 37r J) 



+ . . . , (36) 



where the term e IZ is not relevant for our calculation and has been neglected. Taking the ratio between 
the coefficients of the term gT iz in eqs (f31)|) and (|3"2"|) , we obtain the monodromy to second order: 

M(r ) = -[1 +2COS0V)] (l + -^L= + A2c + M . (37) 

Here, 

_ (l + »e 3 "^ + + (z + e 3 ^- , . 

Al (-1 + e* 3 ^) ; (38) 

A 2c = -(1 -*)£+£_ -£c + _; (39) 

(1 + e i37ri )(d ++ e l7rj / 2 + d__ c^'/ 2 ) - 2d+- - 2d-+ e l3 ^ , . 

= (-1+e^) ' (40) 

The terms A 2c and A 2t ; depend on coefficients and , respectively. Although our expression for 
Ai here is different from that in Ref. [13j by a phase factor, our final result is identical to their. 
Making use of the formula 

f°° v ^7r72r( 1 + 2 ^+ 2v ) 

h{^v) = ^ dzZ-V 2 J„{z)J u {z) = 7^ 3 ^_ 2 , )r( 3+2M _ 2 t )r( 3-2^+2, ^ ( 41 ) 

one can obtain explicitly 

_ (3/ 2 + 3i + i- j 2 ) r 2 (i) r(ifi) r(i±2i) sin[^^] 



48ttI sin(if) 

(3Z 2 + 3/ + 1 - j 2 ) r 2 (±) r(i^i) r(i±2i) sin[«] 



48 tt* sin(^) 



(42) 
(43) 



2 

*-(l-2j)l „:„r^(l+2j)i 



(sz 2 + + 1 - j 2 ) r 2 (i) r(ifi) r(±±^) sin[^M] sin L 



Note that 



= ^ . ■ V4/ v 4 ^__v 4 / * J L 4 £ (44) 

24V27rl sin(^) 



-c++(j -j); c +- = -c-+(j -j)- (45) 



These relation are also obeyed by d^'s, which can be used to reduce our work. With the above results, 
we are ready to find Ai and A 2c in eq (|40]) : 



»(3Z 2 + 3/ + 1 - j 2 ) T 2 (l) T(l^) r(^) cos(f) cos( J7 r) 

A i t^tttttt-t — 7^ si ; ( 4b ) 



6\/27r 3 / 2 [l + 2cos(j7r)] 

(^)r 2 (^ 

11527T 3 



a (3/ 2 + 3/ + 1 - j 2 ) 2 r 4 (i) r 2 (i^) r 2 (i±^) cos( J7 r) 



G 



The double integral 

f°° f Z2 _ 1/2 _ 1/2 
h{H2,V2\lHi v \) = / d2 2 / dziz 2 ^2(22)^2(22)^! (21) Ji/! (21) (48) 

Jo Jo 

can be expressed in terms of the generalized hypergeometric functions, but the general formula is quite 
complicated and not particularly illuminating. Therefore, we will just give the final result explicitly 
for the coefficients d++ and d^ : 

d ^ 2 {3P + 31 + 1-J 2 ) 2 cot(f )r(l) 5 G 4 (l, 11 |, 

++ 576 sin 2 (^1) 

, V5F (3/ 2 + 3/ + 1 - j 2 ) 2 cot(^) cot(j tt) r(i± 2 i) r(ifi) r 2 (i±i) 



288sin(f) 

r ri 1 + 2j 1+j 1 + 3 1 + 2j 2 + j 2 + j 5 + 2j 1 + • 11 

5 4 4' 4 ' 2 ' 2 ' 2 ' 2 ' 2 ' 4 > i+J '- L '' 
, (3Z 2 + 3/ + 1 - ff T\\) r 2 (l^) r 2 (lfl) sin 2 [g^M] singgg] 



1152^2 ir 3 sin 2 (^) 



tt 2 (s/ 2 + 3/ + 1 - j 2 ) 2 r(i) 5 G 4 (i, 1 1 i±i, V; !> 1. t 1 . ¥; 1) 



d-\ = — 



4 ' 2 ' 2 ' 2 ' 2 ' 4 ' 

cin^i'Ll 



1152 sin 3 (^) 

(3/ 2 + 3/ + 1 - j 2 ) 2 cot 2 (^) cot(j tt) r(i^) r(^i) r 2 (i^) 



576 

,1 1 - 2 j 1- j 1-j 1-2,7 . 2- j 2 - j 5-2j _ 
>2' 4 ' 2 ' 2 ' 2 ' 2 ' 2 ' 4 ' 3 

2-LQ/J.1 _,-21 2 r^Il^i^N r 2/l+2is g;T ,2r^(l-2j)l d n 2 r £(1+22)1 



(3l 2 + 3l + l-fy r 4 (i)r 2 (i^)r 2 (l±^) sin^^^l] sin 



4 



2304 tt 3 sin 2 (^) 



(50) 



Here, we have used the regularized generalized hypergeometric function 56*4(01, a%, a$, CI4, 05; 61, 62, &3, &4! 2) 
so that the pole structure of each term in these expressions are more explicit. It is related to the usual 
generalized hypergeometric function by 

n , , , , , s 5-^4(01,02, 03, a 4 ,,a 5 ;b 1 ,b 2 ,b 3 ,b 4 ; z) 
5 G 4 (01,02,03,04, a 5 ; b u b 2l b 3l b 4 ; z) = — . 51 

r(6i)r(&2)r(6 3 )r(64) 

The other two coefficients can be obtained by relations analogous to those in eq (|4"S")l 

d— = -d++(j -f -i); d_+ = -d+_(i^-i). (52) 

On the face of it, each of the d^'s has a third order pole coming from terms involving the generalized 
hypergeometric function when j is an even integer. On closer look, we see there are some cancelation 
among the divergences and in the end all they have are just simple poles in such limit similar to the 

c^j/'s. Another possible divergence arises in d when j = 1, which will again be canceled when we 

calculate the monodromy. 
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It is now straightforward to obtain &2d by making use of the following two identities 



-4^/ 2 cos2( ^) eot( J7 r)r(i±^)r(i±^)r 2 (i±^) 

r (1 1 + 2 -7 1+1 1+1 1 + 2 -? 2 + J 2 + J 5 + 2 J , , ■ u 
5U4 4' 4 ' 2 ' 2 ' 2 ' 2 ' 2 ' 4 

+47r V2 cos2( ^ } cot(j7r)r( i_^)r(i^)r 2 (i y ^) 



f J J_ J_ J . J -J_ -J i _ ■ 1 \ 

^ 2 ' 4 ' 2 ' 2 ' 2 ' 2 ' 2 ' 4 ' 

-r^)r^)r^) ^ sin[^±M] = o ; (53) 
-4.^ cos( |) r ( i±^)r ( i±^V(i±i) 
r ri 1 + 2j l±i 1±2 1 + 2 -? 2 + J 2 + J j+jj . 

5 4l 2' 4 ' 2 ' 2 ' 2 ' 2 ' 2 ' 4 

cos(f )r(i^)r(i^)r 2 (i^) 

r rl LljZ L^j !-J i-2j 2- j 2- j 5-2j 
5 4l 2' 4 ' 2 ' 2 ' 2 ' 2 ' 2 ' 4 ' h h 
,o5 r A n m 1+J 1-J.5 1 2+j 2-j 

-r 4 (i)r 2 (i^)r 2 (l±^) cos (f )[i-cos(i 7r )] = o. (54) 

Eventually, we achieve the following nice result 

a _ (3^ 2 + 2,1 + 1- j 2 ) 2 T\\) r 2 (l^) r 2 (^) cosQ-tt) 
2d 1152^ 3 [l + 2cos(j7r)] ' 1 ' 

where all divergences have been canceled out. 

Together with the result from eq (|47[) . the asymptotic form of quasinormal frequencies of a four 
dimensional Schwarzschild black hole is found to be 

47iw = (2n + l)iri + ln[l + 2cos(j7r)] 

^(3^ 2 + 31 + 1- j 2 ) T 2 (l) r(l±2i) cos(f ) cos(jTr) 

%\f^Xiv > 1 2 [1 + 2cos(j7r)]V^^ 

(3? 2 + 3/ + 1 - j 2 ) 2 r 4 (i) r 2 ( r 2 ( ifi) CO s 2 (^) 

576tt 3 [1 + 2 cos(j7r)] 2 (-wr ) 
The physically interested cases are 

co n i - z (z 2 + z - i)r 4 (i/4) i (/ 2 + z- i) 2 r 8 (i/4) 



+ O[(-^r )- 3 / 2 ]. (56) 



— w (2rc + l)7ri + ln3 , . 
T H 1 7 yft 18V27T 3 / 2 n 2592^ 3 

W "W9 :_ui o i l-^ (^ 2 + ^ + l/3)r 4 (l/4) » (i 2 + i + l/3)¥(l/4) 

fa zi + 1 hrz + m3 H = = — ; - . 

T H K ' ^ 6\/2V/ 2 n 288tt 3 

2mti 



u n „ . i2n(l 2 + 1) 2 



for j = 2; 


(57) 


for j = 0; 


(58) 


for j = l. 


(59) 



T H n 

A few comments are in order. First, all the second order corrections are purely imaginary. In particular, 
when j = 2 (gravitational perturbation) the numerical coefficients of the i/n term (after divided by 
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Att) are 0.739,3.58,49.7 for I — 2,3,6, respectively. They are in good agreement with the known 
numerical studies [3]. As for the real part, our result predicts vanishing correction. For j = 2, this is 
again consistent with the numerical results in Ref. [I] for I — 2,3. For I = 6 the numerical result is 
0.263, which seems to be contradictory to ours. However, the numerical value for I = 6 has opposite 
sign relative to those of I = 2, 3. This is peculiar, since in all other cases a given type of corrections 
are always of the same sign irrespective of the specific value of angular momentum. Therefore, we 
believe more study is needed to clarify whether there is really a discrepancy. As for the j = 1 case, 
the numerical study in Ref. [T7] suggests the leading correction is of the form —^j^. However, this 
does not necessarily mean the two results are inconsistent. In fact, one can only extract the behavior 
of the leading correction to the real part from their Fig. 2 and further numerical study is needed to 
confirm or refute our prediction. 

3 Conclusion 

In sum, we have calculated to second order the correction to the asymptotic form of quasinormal 
frequencies for Schwarzschild black holes in four dimensions. Most of our results are consistent with 
the numerical ones when available. In cases where there seem to be contradiction, we think further 
numerical studies are needed to clarify the situation. It would also be helpful if more detailed numerical 
studies can be carried out for the j = case so that more thorough comparisons are possible. It would 
be interesting to generalize the method to other spacetime backgrounds [TB] . Extension to higher order 
is also desirable. It might enable us to find a quantitative prediction for the "algebraically special" 
frequencies in Schwarzschild black holes, where the quasinormal frequency is purely imaginary and it 
increases with the fourth power of I [19l [4] . 
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